Charge Kondo anomalies in PbTe doped with Tl impurities 
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We investigate the properties of PbTe doped with a small concentration x of Tl impurities acting 
as acceptors and described by Anderson impurities with negative onsite correlation energy. We use 
the numerical renormalization group method to show that the resulting charge Kondo effect natu- 
rally accounts for the unusual low temperature and doping dependence of normal state properties, 
including the self-compensation effect in the carrier density and the nonmagnetic Kondo anomaly 
in the resistivity. These are found to be in good qualitative agreement with experiment. Our results 
for the Tl s-electron spectral function provide a new interpretation of point contact data. 
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tures: while for x < x c , the residual resistivity, po, is very 
small and almost constant as a function of x, for x > x c , 
po increases linearly with x Q. For x < x c , the resistiv- 
ity, p{T), exhibits a positive slope at low temperature, 
while for x > x c , the slope is negative and a Kondo-like 
contribution, pi mp (T), is observed for T < 10K 0, Q. 
The origin of this anomaly is not due to magnetic impu- 
rities, since the susceptibility is diamagnetic @. 

In this Letter, we show that the unusual concentra- 
tion and low temperature dependence of a number of 
normal state properties of Pbi-^Tl^Te can be naturally 
explained within a picture of dilute Tl impurities acting 
as negative U centers in the PbTe host. Such a dilute im- 
purity description is suggested by the observed concen- 
tration dependence of many properties, as summarized 
above, e.g., the linear dependence on x of the residual 
resistivity for x > x c . We model the Tl impurities by 
the negative U Anderson model [l3| and solve this for 
the static and dynamic properties by using the numerical 
renormalization group (NRG) method [15j. Our results, 
with comparisons to experimental data, and previous ex- 
perimental and theoretical work on Pbi-^Tl^Te [tj.l7l.fl6l]. 
provide strong evidence that the CK effect [3| is realized 
in this system. 

The idea that doping PbTe with Tl induces resonant 
states in the valence band of PbTe was conjectured early 
on (seeRef.Q). These states are mainly of Tl s-character 
and lie close to the top of the valence band, as re- 
cently shown by density functional theory calculations 
17 1. Resonant states alone, however, and generaliza- 
tions of this to an impurity band of resonant states, can- 
not explain properties such as the superconductivity of 
Pbi _ x Tl x Te. For this, a coupling of Tl ions to the lattice 
[lq . or a static mixed valence model [l^| have been pro- 
posed. In the latter, Tl impurities, known to be valence 
skippers in compounds, are assumed to dissociate into 
Transport mea- energetically close Tl 1+ (6s 2 p 1 ) and Tl 3+ (6s°p 3 ) ions, 
surements also show anomalous behavior at low tempera- while the Tl 2+ (6s 1 p 2 ) configuration lies higher in energy 



Introduction. — PbTe is a narrow gap semiconductor 
with a band gap of 190 meV at zero temperature [l|. 
Upon doping with Tl impurities 0, Q a number of strik- 
ing anomalies appear in the low temperature properties. 
These anomalies, described below, are currently not well 
understood. They are believed to be due to the spe- 
cial nature of Tl impurities in PbTe, which have been 
proposed to act as negative U centers, where U < is 
an attractive onsite Coulomb energy. In compounds, the 
outer 6 s shell of Tl can be empty, singly occupied or dou- 
bly occupied. A negative Coulomb correlation results in 
the nonmagnetic empty and doubly occupied states being 
lower in energy than the magnetic singly occupied state. 
Hybridization of the s-states with the valence band states 
of PbTe can then result, via dynamic valence fluctua- 
tions, in a (nonmagnetic) charge Kondo (CK) analogue 
0, Hj] of the conventional spin Kondo effect This 
nonmagnetic CK effect has been argued to explain many 
of the observed properties of Pbi-^Tl^Te [6 -flo| which 
would then constitute the first physical realization of this 
effect. 

A remarkable feature of Pbi-^Tl^Te is that the low 
temperature properties depend sensitively on the Tl con- 
centration, with a qualitatively different behavior below 
and above a critical concentration x c — 0.3 % Tl. For 
example, while Pbi_ x Tl x Te remains metalliclike down 
to the lowest temperatures for x < x c , for x > x c it 
becomes superconducting with a transition temperature 
T c (x) increasing linearly with x and reaching 1.5 K at 
x = 1.5 % Tl 0, [ll[. This is surprisingly high given the 
low carrier density of less than 10 20 holes/cm 3 . In ad- 
dition, measurements of the Hall number = l/R^e 
[E EH indicate that the number of holes grows linearly 
with x for x < x c , whereas for x > x c the number of 
holes remains almost constant: the system exhibits "self- 
compensation" and the chemical pot ential is pinned to 
a value p 
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[20|. In the highly polarizable PbTe host, this can re- 
sult in negative on-site U and provides a mechanism for 
superconductivity and a qualitative explanation for the 
observed self-compensation, chemical potential pinning 
and the diamagnetic behavior of Pbi^TLjTe [19(. For 
a more quantitative explanation of the observed anoma- 
lies, and in order to explain the Kondo anomalies in the 
resistivity, a more realistic model is needed, which in- 
cludes dynamic fluctuations between the Tl 1+ and Tl 3+ 
valence states. This motivates our use of the negative 
U Anderson model [14(, as formulated for Pbi-^Tl^Te 



in Ref. M and discussed as a model for this system in 
Refs. la-llOi lid l2lL For completeness, we mention also the 
valence band model for PbTe, in which the main effect of 
Tl doping is assumed to be a rigid shift of the chemical 
potential into the valence band. This model may be rele- 
vant for the transport properties of PbTe doped with Tl 
impurities 22j at high temperature (T > 300i^) , where 
the charge Kondo effect is suppressed. It fails, however, 
to describe the low temperature anomalies that we are 
addressing in this Letter, e.g. the Kondo upturn in the 
resistivity at T < 10 K for x > x c . 

Model and calculations. — We consider n Tl impuri- 
ties in a PbTe crystal with N Pb sites described by the 

H imp + H hyb . The first term, 



Hamiltonian H = H. 



band 



Hband = Z)ko-( e k _ Ue)cl a c kcri describes the valence p- 
band of PbTe, where u e is the (electron) chemical poten- 
tial and c^ a creates an electron with energy ek- The sec- 
ond term, H imp = (e -fx e ) J2i=ia n isa+U Y!i=i n ls \n lsil 
describes the Tl impurities, where ni sa = s\ a s ia is the 
number operator for a Tl s-electron at site i with spin a 
and energy eo, and U is the (negative) correlation energy. 
The last term, H hyb = Yh=i J2ka Voic^s^ + h.c), mod- 
els the hybridization of Tl s-states with the valence band 
p-states and Vq is the matrix element for the s-p inter- 
action. Its strength is characterized by the hybridization 
function A(w) = ttVq J2k — £ kcr) = nV Q 2 Af(u}), where 
we retain the full energy dependence of the p-band den- 
sity of states Af(w) = Ek 6 ( UJ ~ £ k) of PbTe [H]. 

The chemical potential u e determines n e = 
^Ek<r( c L c k<r) and n s = ^^(fli,,), the aver- 
age number of p and s electrons per site. We denote by 
x = n/N the concentration of Tl impurities. Since Tl 
acts as an acceptor, the ground state corresponds to the 
Tl 1+ (n s = 2) configuration and the Tl 3+ (n s = 0) con- 
figuration is split-off from the ground state by the energy 
5 = E(Tl a+ ) - E(Tl 1+ ) > 0. A concentration x of Tl im- 
purities accommodates x(n s — 1) electrons (per Tl site), 
where the number of accepted electrons in the 6s level 
of Tl is measured relative to the neutral Tl 2+ (s 1 ) con- 
figuration having n 8 = 1. These electrons are removed 
from the valence band leaving behind no = 1 — n e holes. 
Thus, charge neutrality implies no = x{n s — 1) Q, which 
for a given x and temperature T has to be satisfied by 
adjusting the chemical potential u e . Here, we neglect in- 



tcrimpurity interactions and solve H for a collection of 
sing le independent negative-C/ centers by using the NRG 
[15J . For each x and each T we satisfy the above equation 
by sclf-consistently determining the chemical potential « e 
(or equivalently the hole chemical potential /i, which we 
henceforth use). 

The electrical resistivity of electrons scattering from a 
dilute concentration x of Tl impurities is obtained from 
the usual expression, pi m p{T) = l/e 2 Ln, where L\\ is 
the static limit of the current-current correlation func- 
tion. In the absence of nonresonant scattering the vertex 
corrections vanish and the relevant transport integral L\\ 
can be written as 1231. 
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dm 



M(u)t(u,T), (1) 



where 00 = ( v k F ) is the velocity factor v\ averaged over 
the Fermi surface, f(ui) = 1/[1 + cxp(uj /ksT)] is the 
Fermi function, t{ui,T) is the conduction-electron trans- 
port timed! r(w,T)~ 1 = 2n T \Vg A{uj ,T) and A(ui,T) = 
— —lm G(u> + i0 + ) is the spectral function of s-electrons. 
The number of Tl impurities nxi per cm 3 is related to x 
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in percent by nxi = 1.48 x 10 

Choice of model parameters. — We use u* w 225 meV, 
close to the value obtained from tunneling experiments 
[lil . 13 1 . Other parameters, such as A = A(u*), re- 
quired to fix the hybridization function A(cj), and U 
are largely unknown. Our interpretation of the tun- 
neling spectra (see discussion of Fig. [2] below), suggests 
U w — 30meV. The measured Kondo-like resistivity for 
x > x c at low temperatures requires that \U\ 3> Ao- We 
take U/Aq = —11 with Ao = 2.7 meV to yield a Kondo 
temperature Tjc w 1.23 K below the highest T c « 1.5 K 
at x = 1.5% 3> x c , where Tk is defined via the im- 
purity resistivity p lmp (T = Tk) = 0.5pi mp (T = 0) for 
x x c - The overall qualitative aspects of our results 
remained the same for values of U in the range lOmeV 
< \U\ < 220meV and \U\/A > 1. 

Qualitative considerations. — It is instructive to first 
make some qualitative remarks, starting from the atomic 
limit Vq = 0. For x = the chemical potential lies 
in the gap between the valence and conduction bands. 
For finite but very small x each Tl impurity accepts one 
electron, i.e. n s sa 2 and no = x(n s — 1) x grows lin- 
early with x. At the same time, the chemical potential 
shifts downwards into the valence band u < E v , where 
E v denotes the top of the valence band. This implies that 
the splitting S(u) = — (2(eo — u) + U) between donor and 
acceptor configurations decreases. Eventually, at a criti- 
cal concentration x = x* , the chemical potential reaches 
u = /i* = £ + U/2 where 8(u) = 2(u — (j,*) = and 
the system is in a (static) mixed valence state where the 
Tl 1+ and Tl 3+ configurations are degenerate. In this sit- 
uation n s = 1, and further doping cannot increase the 
hole carrier density beyond the value no (//*), i.e. one has 
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FIG. 1: (Color online) Hole carrier density no versus Tl dop- 
ing a: in % for T = 1.8 K and T = 77 K. Filled circles: exper- 
imental data at T = 77 K filled squares: experimental 
data at T = 1.8 K dot-dashed line: expected no for one 
hole per Tl. Inset (a): hole chemical potential li versus x at 
T = 1.8K and T = 77K [and hence 5(fx) = 2(/i - u*))]. 



self-compensation with a pinning of the chemical poten- 
tial to LI* 0]. 

For finite Vq , quantum fluctuations between the degen- 
erate states Tl 1+ and Tl 3+ at fx = fx* become important 
and lead to a CK effect. This significantly affects all 
static and dynamic properties and needs to be taken into 
account in describing the experiments. It is also impor- 
tant for fi > fi* , since a finite charge splitting 5(ja) > 
in the negative-?/ Anderson model is similar to a Zee- 
man splitting in the conventional spin Kondo effect [25| . 
The latter is known to drastically influence all proper- 
ties Thus, for the whole range of concentrations x, 
one expects fluctuations to play an important role in the 
properties of Pbi-^Tl^Tc. 

Numerical results. — Figure [T] shows rio(x) versus x at 
T = 1.8 K and at T = 77 K and a comparison with exper- 
imental data on Hall number measurements [H, 12 1. At 
low dopings, no is linear in x both in theory and in exper- 
iment, as expected for Tl impurities acting as acceptors 
(dot-dashed curve in Fig. [1]). However, the efficiency, 
no(x)/riTi(a;), of Tl dopants in supplying holes at low x 
is only around 65 % in the data of Ref. [12J as opposed to 
100 % in our model calculations and in the data of Ref. H|. 
At higher dopings no (a;) saturates rapidly with increas- 
ing x for T = 1.8 K and more slowly at higher tempera- 
tures. The theoretical crossover from linear to saturated 
behavior occurs at x = x* ps 0.5 %, larger than the value 
x c re 0.3 % for the onset of superconductivity. The the- 
oretical saturation density rip re 0.7 x 10 20 /cm 3 is close 
to the experimental value [12j. The self-compensation 
effect for x 3> x* re 0.5% is a characteristic signature 
of the CK state: on entering this state the Tl ions fluc- 
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FIG. 2: (Color online) Tl spectral function A(E, T = 0) ver- 
sus E — Ef for a range of Tl dopings x. Inset (a): region near 
E = Ef showing the CK resonance. Inset (b): PbTe valence 
band densit y o f states N(E) versus E — Ef for x as in the 
main panel [22| (legend: hole densities no for each x). 



tuate between Tl 1+ and Tl 3+ so the average valence is 
Tl 2+ , which corresponds to no additional electrons being 
accepted or donated. The CK state may be destroyed 
by lifting the degeneracy of the pseudospin states, i.e. 
in the language of the spin Kondo effect, by applying a 
"magnetic field" . A "magnetic field" in the CK effect cor- 
responds to doping or shifting fx. This has been achieved 
by counterdoping with In ions (lo| . which act as donors. 
The Kondo anomalies, e.g. in the resistivity, were ob- 
served to vanish, providing further support to the CK 
picture. The self-compensation effect can be seen also in 
the pinning of the hole chemical potential \i for x > x* , 
shown in Fig. QJi. For x < x* , \i grows nonlinearly with 
x and rapidly approaches the value fi* = 225 meV for 
x > x*, both at T = 1.8K and at T = 77K. 

Figure [2] shows the Tl s-electron spectral function 
A(E, T = 0) at zero temperature and different dopings 
x. For small doping, x < x* , the hole chemical potential 
Li lies above li* within the shallow part of the valence 
band density of states N(E) (see Fig. [2Jd), consequently 
the splitting 8{ll) = 2(fj, — fj,*), see Fig. QJi, is large. Such 
a large splitting acts like a large Zeeman splitting in the 
conventional positive U Anderson model and polarizes 
the spectral function so that its weight lies mostly below 
the Fermi level Ep For x > x* , li approaches li* and 
a CK effect develops. The spectral function develops a 
sharp asymmetric Kondo resonance close to, but below 
Ef (see Fig. [2k,) and an upper Hubbard satellite peak 
appears above Ep- Early tunneling experiments for low 
dopings x < 0.3% showed only one resonant level be- 
low Ep 13j, whereas more recent tunneling experiments 

a nar- 
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for x > 0.6 % show two "quasilocal" peaks 
row one of width 6 meV close to Ep and a broader one 
of width 12meV at a nearly constant energy 13-15 meV 
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FIG. 3: (Color online) Impurity resistivity pi mp (T) / pu versus 
temperature T and a range of Tl concentrations x (in percent, 
legend: column 1). For x > x* , pj mv is well described by the 
usual spin Kondo resistivity [23] (shown as open squares for 
x = 1.3%), but with effective Kondo scales T^ff (legend: 
column 2). p u = 2titi / (e 2 TvfiNpao) is the residual resistivity 
for unitary scatterers with Nf = Af(p*). The vertical arrow: 
7k = 1.23 K. Inset (a): p = xpi mp (T — 0)/p u versus x (filled 
circles). Dashed line: pa for static unitary scatterers in place 
of Tl. 



below this. Interpreting these as the Kondo resonance 
and the lower Hubbard satellite peaks in Fig. [2] yields 
\U\/2 « 15meV and hence the value U « — 30meV used 
in our calculations. Since the CK resonance is temper- 
ature dependent the above interpretation could be 
tested with temperature dependent studies of tunneling 
or high resolution photoemission spectra. 

The temperature and doping dependence of the impu- 
rity resistivity, pi mp , is shown in Fig. [3] For x < x* , Tl 
impurities act as acceptors with a well-defined valence 
state (Tl 1+ ). They therefore act as weak potential scat- 
terers and consequently the resistivity is much below the 
unitary value, as seen in Fig- El For x > x* , dynamic fluc- 
tuations between the nearly degenerate Tl 1+ and Tl 3+ 
states leads to the CK effect and pi m p approaches the 
resistivity for unitary scatterers at T = 0. For x > x* , 
Pimp is well described by the spin Kondo resistivity [23J 
(open squares, Fig. [3]) with a logarithmic form around 
T w T^ ff , where T e J f is an effective Kondo scale, and 
T 2 Fermi liquid corrections at low T -C T^f^, in qualita- 
tive agreement with experiment @ . The effective Kondo 
scale Tftf is a function of the charge splitting S(p) and 
Tr-, and approaches the true Kondo scale Tk only asymp- 
totically for x » x* (see legend to Fig.[3j. Finally, Fig.[3tt 
shows that the impurity residual resistivity is significant 
only when the CK effect is operative, i.e. for x > x* , in 
qualitative agreement with experiment @, Q . 

Conclusions. — In summary, we investigated the nor- 
mal state properties of Pbi-^Tl^Te within a model of 



Tl impurities acting as negative U centers. Our NRG 
calculations support the suggestion that the CK effect is 
realized in Pbi-xTl^Te @, 13|- They explain a number 
of low temperature anomalies of Pbi-^Tl^Te, including 
the qualitatively different behavior below and above the 
critical concentration x* , where x* sw 0.5% is close to 
x c m 0.3 % for the onset of superconductivity. At x = x* , 
two nonmagnetic valence states of Tl become almost de- 
generate and the ensuing pseudospin CK effect results in 
a Kondo anomaly in the resistivity for x > x* and a resid- 
ual resistivity approximately linear in x. Our results for 
these quantities and the carrier density nntsc) are in good 
qualitative agreement with experiments SSSO- For 
the Tl s-electron spectral function, we predict that one 
peak should be present far below Ep for x < x* and that 
a second temperature dependent Kondo resonance peak 
develops close to, but below Ep, on increasing x above 
x* . This provides a new interpretation of measured tun- 
neling spectra [T^j, which could be tested by temperature 
dependent studies of tunneling or photoemission spectra. 
In the future, it would be interesting to extend this work 
to include the effects of disorder within an Anderson- 
Hubbard model description. 
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